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Abstract 

Two computer codes relevant for the description of deep inelastic scattering off polarized 
targets are discussed. The code uela deals with radiative corrections to elastic fie scattering, 
one method applied for muon beam polarimetry. The code HECTOR allows to calculate both 
the radiative corrections for unpolarized and polarized deep inelastic scattering, including 
higher order QED corrections. 



1 Introduction 



The exact knowledge of QED, QCD, and electroweak (EW) radiative corrections (RC) to the 
deep inelastic scattering (DIS) processes is necessary for a precise determination of the nucleon 
structure functions. The present and forthcoming high statistics measurements of polarized 
structure functions in the SLAC experiments, by HERMES, and later by COMPASS require the 
knowledge of the RC to the DIS polarized cross-sections at the percent level. 

Several codes based on different approaches for the calculation of the RC to DIS experiments, 
mainly for non-polarized DIS, were developped and thoroughly compared in the past, cf. [|T|. 
Later on the radiative corrections for a vast amount of experimentally relevant sets of kinematic 
variables were calculated ||, including also semi-inclusive situations as the RC's in the case of 
tagged photons ||. Furthermore the radiative corrections to elastic \i-e scattering, a process to 
monitor (polarized) muon beams, were calculated M. The corresponding codes are : 

• HECTOR 1.00, (1994-1995) §, by the Dubna-Zeuthen Group. It calculates QED, QCD 
and EW corrections for variety of measuremets for unpolarized DIS. 

• fxela 1.00, (March 1996) J|], calculates 0(a) QED correction for polarized /ie elastic 
scattering. 
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HECTOR 1.11, (1996) extends HECTOR 1 . 00 including the radiative corrections for polarized 
DIS ||, and for DIS with tagged photons 0. The beta- version of the code is available 



from pittp : // www . if h . de/ 



2 The Program fiela 

Muon beams may be monitored using the processes of fi decay and fie scattering in case of 
atomic targets. Both processes were used by the SMC experiment. Similar techniques will be 
used by the COMPASS experiment. For the cross section measurement the radiative corrections 
to these processes have to be known at high precision. For this purpose a renewed calculation 
of the radiative corrections to cr(/xe — > fie) was performed ||. 

The differential cross-section of polarized elastic \xe scattering in the Born approximation 
reads, cf. M, 



, BORN rj 

da Ziia 



dy m e E l 



2^ (1 - yP e P,) + ± (1 - P e P,) 



(1) 



where y = y {t = \- E'^/E^ = E'jE^ = y e , Y = (1 + m^/2/E^) 1 = y max , m M , m e - muon and 
electron masses, E^, E„, E e the energies of the incoming and outgoing muon, and outgoing 
electron respectively, in the laboratory frame. P M and P e denote the longitudinal polarizations 
of muon beam and electron target. At Born level y^ and y e agree. However, both quantities are 
different under inclusion of radiative corrections due to bremsstrahlung. The correction factors 
may be rather different depending on which variables (y^ or y e ) are used. 

In the SMC analysis the ^-distribution was used to measure the electron spin-flip asymmetry 
A e *v. Since previous calculations, || |9|, referred to y e , and only ref. || took polarizations into 
account, a new calculation was performed, including the complete 0(a) QED correction for the 
^-distribution, longitudinal polarizations for both leptons, the //-mass effects, and neglecting 
m e wherever possible. Furthermore the present calculation allows for cuts on the electron re- 
coil energy (35GeV), the energy balance (40GeV), and angular cuts for both outgoing leptons 
(1 mrad). The default values are given in parentheses. 

Up to order 0(a 3 ), 14 Feynman graphs contribute to the cross-section for \i-e scattering, 
which may be subdivided into 12 = 2 x 6 pieces, which are separately gauge invariant 

j Q ED 2 6 i i 

r,n ^2^2 — (2) 



% i=i k =i % 



One may express (E3) also as 



r QED / , url pol j pol \ 

^ = yf d(J k + da% \ _ 

dy^ ^ V d Vv M d Vi* J ' 

The indices I and k in the combinations Ik have the meaning 

1=1 — unpolarized contribution, I = unpol; 

2 — polarized contribution (terms with P e P M in g?<t born ), I = pol. 
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k = 1 — Born cross-section, k = b; 

2 — RC for the muonic current: vertex + bremsstrahlung, k = fifi; 

3 — amm contribution from muonic current, k = amm; 

4 — RC for the electronic current: vertex + bremsstrahlung, k = ee; 

5 — fie interference: two-photon exchange + 

muon-electron bremsstrahlung interference, k = fie; 

6 — vacuum polarization correction, running a, k = vp. 

The FORTRAN code for the scattering cross section (|2|) fiela was used in a recent analysis of the 
SMC collaboration. 

The RC, 8y , to the asymmetry A^ u shown in figures 1 and 2 is defined as 

< = (%). where A, e = (4) 

The results may be summarized as follows. The 0(a) QED RC to polarized elastic fie scattering 
were calculated for the first time using the variable y^. A rather general FORTRAN code fiela for 
this process was created allowing for the inclusion of kinematic cuts. Since under the conditions of 
the SMC experiment the corrections turn out to be small our calculation justifies their neglection. 



3 Program HECTOR 

3.1 Different approaches to RC for DIS 

The radiative corrections to deep inelastic scattering are treated using two basic approaches. 
One possibility consists in generating events on the basis of matrix elements including the RC's. 
This approach is suited for detector simulations, but requests a very hughe number of events 
to obtain the corrections at a high precision. Alternatively, semi-analytic codes allow a fast 
and very precise evaluation, even including a series of basic cuts and flexible adjustment to 
specific phase space requirements, which may be caused by the way kinematic variables are 
experimentally measured, cf. 0, |5|. Recently, a third approach, the so-called deterministic 
approach, was followed, cf. |I(J. It treats the RC's completely exclusively combining features of 
fast computing with the possibility to apply any cuts. Some elements of this approach were used 
in fiela and in the branch of HECTOR 1.11, in which DIS with tagged photons is calculated. 

Concerning the theoretical treatment three approaches are in use to calculate the radiative 
corrections: 1) the model-independent approach (MI); 2) the leading-log approximation (LLA); 
and 3) an approach based on the quark-parton model (QPM) in evaluating the radiative correc- 
tions to the scattering cross-section. 

In the model-independent approach the QED corrections are only evaluated for the leptonic 
tensor. Strictly it applies only for neutral current processes. The hadronic tensor can be dealt 
with in its most general form on the Lorentz- level. Both lepton-hadron corrections as well as 
pure hadronic corrections are neglected. This is justified in a series of cases in which these 
corrections turn out to be very small. The leading logarithmic approximation is one of the 
semi-analytic treatments in which the different collinear singularities of 0((aln(Q 2 /mf)) n ) are 
evaluated and other corrections are neglected. The QPM-approach deals with the full set of 
diagrams on the quark level. Within this method, any corrections (lepton-hadron interference, 
EW) can be included. However, it has limited precision too, now due to use of QPM-model itself. 
Details on the realization of these approaches within the code HECTOR are given in ref. || |TT|| . 
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3.2 0(a) QED Corrections for Polarized Deep Inelastic Scattering 

To introduce basic notation, we show the Born diagram 




p(p,M) X(p',M h ) 

and the Born cross-section, which is presented as the product of the leptonic and hadronic tensor 

2na 2 



da 



Born 



with 

q = ki- k 2 , 
and the Bjorken scaling variables 



V 



Q 4 

Q 2 

p.q 



-y 



dxdy, 
S = 2(p.h) 



x 



Q 2 



p.kt Sy 
For the hadronic tensor, we use the representation of ref. [|12|] 

p^2 7 r) 6 ^/(p'|^b)(p|Xb> 4 (Eft-p')n^ 



(5) 
(6) 

(7) 



W, 



-9uu + 



q- 

q x s a 



+ieau\a- Gi(x, Q ) iie^xiT 

p.q 



p.q 

q x (p.qs a — s.qp a 



q x p u 



(p.qf 



2p.q 

G2(x,Q 2 



^s(x,Q 2 



+ 



P^u + S^p v PnP u 

s.q-^- 

2 p.q 



—G 3 (x,Q 2 

p.q 



+s.q 7 J ^ 2 -Gi{x, Q ) + -g^ + — y s (x, Q 



(p.q) 



Q ) V-1 



where 



Pfi — Pa ^ 2 Qn, 



s fj, — S H 2 



S4 

q 2 



and s is the four vector of nucleon polarization, which is given by s = A p M(0, n) in the nucleon 
rest frame. 
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The combined structure functions in eq. (§) 



^i, 2 (x,Q 2 ) 
£3,4,5(2, Q 2 ) 



Q 2 e F?J( Xl Q 2 ) + 2\Q e \ ( Vl - Pe \iai)x(Q 2 )Fl z 2 (x,Q z 

}2\ 



+ + a 2 - 2p e A M a,) x\Q 2 )F?,i& Q 2 
2|Q e |(p e a i -V0x(Q 2 )F 3 ^(x,Q 2 ), 

2p e ^a ; - A; (v? + a 2 )] X 2 (Q 2 )Fi z (x, Q 2 



+ 



-Q£>*9l[l(x, Q 2 ) + 2|Q e | (p e a, - A M ) x(Q 2 )<?Z.I Q 



+ 



2p e ^ - A ; (v? + a 2 )] x 2 (g 2 )^ 2 z (x, Q 2 ), 

2|Q e | (ui - PeAjaj) xCQ 2 )^^, <5 2 )> 

+ (t; 2 + a 2 - 2p e W,) xWK^, Q 2 ), (9) 

are expressed via the hadronic structure functions, the Z-boson-lepton couplings Uj, Oj, and the 
ratio of the propagators for the photon and Z-boson 

Q 2 



x(Q 2 



G,M 2 



T2 ■ 



(10) 



V2 87ra Q 2 + M\ 

Furthermore we use the parameter p e for which p e = 1 for a scattered lepton and p e = — 1 for 
a scattered antilepton. The hadronic structure functions can be expressed in terms of parton 
densities accounting for the twist- 2 contributions only, see |]12|| . Here, a series of relations between 
the different structure functions are used in leading order QCD. 
The DIS cross-section on the Born-level 

d CTBorn 



dxdy dxdy 
contains two contributions, the unpolarized part 



,2 unpol ,2 pol 
u °Born _|_ u °Born 



dxdy 



j unpol 
aa Born 



dxdy 

with the kinematic functions 

sUviQ 2 ) 
S%(y,Q 2 ) 

S%(y,Q 2 ) 



2 -^sj: sr(x,y)H^Q 2 
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(12) 



2a;?/ z 



5 



x 



and the polarized part 



Born 



dxdy 



Q 4 



X p N f p Sj2 S p i (x,y)G l (x,Q 2 



(13) 



(14) 



i=l 



Here, S p gi (x,y) are functions, similar to 
abbrevations 



), and may be found in ||. Furthermore we used the 



f 



a 



1. 



n 



\ki\ 



cos 



4M 2 x ( M 2 xy 
*K\~Sy~\~ V ~~S~ \ 

\ N n±, with ki.n± = 0. 



COS( 



2vr y 



sin 0. 



2- 



(15) 
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The 0(a) DIS cross-section reads 



j2 j2 j2 rfi* un P°l A2 P°l 

« 0"QED,1 _ Ot ^ d 0"Born « 0"Brems _ « ^QED,! " ^QED,! /^g\ 

dxidyi 7r dxidyi dxidyi dxidyi dxidyi 

All partial cross-sections have a form similar to the Born cross-section and are expressed in terms 
of kinematic functions and combinations of structure functions. In the 0(a) approximation the 
measured cross-section, <7 ra d, is define as 

A2 j2 j2 j2 unpol t2 pol 

« 0"rad _ « C"Born « CTQED.l _ « <?Vad _|_ " a rad 

dxidyi dxidyi dxidyi dxidyi dxidyi 

In the four following figures we illustrate the RC-factor 

S= ***__ L (18) 

« CBorn 

The radiative corrections calculated for leptonic variables grow towards high y and smaller values 
of x. The figures compare the results obtained in LLA, accounting for initial (i) and final state 
(/) radiation, as well as the Compton contribution (c2) with the result of the complete calculation 
of the leptonic corrections. In most of the phase space the LLA correction provides an excellent 
description, except of extreme kinematic ranges. 

A comparison of the radiative corrections for polarized deep inelastic scattering between the 



codes HECTOR and POLRAD [17| was carried out. It had to be performed under simplified conditions 



due to the restrictions of POLRAD. Corresponding results may be found in |lT| , 131 14] 



3.3 Conclusions 

For the evaluation of the QED radiative corrections to deep inelastic scattering of polarized 
targets two codes HECTOR and POLRAD exist. The code HECTOR allows a completely general study 
of the radiative corrections in the model independent approach in 0(a) for neutral current reac- 
tions including Z-boson exchange. Furthermore, the LLA corrections are available in 1st and 2nd 
order, including soft-photon resummation and for charged current reactions. POLRAD contains 
a branch which may be used for some semi-inclusive DIS processes. The initial state radia- 
tive corrections (to 2nd order in LLA + soft photon exponentiation) to these (and many more 
processes) can be calculated in detail with the code HECTOR, if the corresponding user-supplied 
routine USRBRN is used together with this package. This applies both for neutral and charged 
current processes as well as a large variety of different measurements of kinematic variables. 
Aside the leptonic corrections, which were studied in detail already, further investigations may 
concern QED corrections to the hadronic tensor as well as the interference terms. 
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Figure 1: The QED radiative corrections to asymmetry without experimental cuts 
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Figure 2: The QED radiative corrections to asymmetry with experimental cuts. 
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Figure 3: A comparison of complete and LLA RC's in the kinematic regime of HERMES for neutral 
current longitudinally polarized DIS in leptonic variables. The polarized parton densities |15[ are 
used. The structure function g 2 is calculated using the Wandzura-Wilczek relation. c2 stands for the 
Compton contribution, see [0 for details. 
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Figure 4: The same as in fig. ^|, but for energies in the range of the SMC-experiment. 
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Figure 5: The same as in fig. |] for x = 10" 
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Figure 6: A comparison of complete and LLA RC's at HERA collider kinematic regime for neutral 
current deep inelastic scattering off a longitudinally polarized target measuring the kinematic variables 
at the leptonic vertex. 
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